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On the distribution of the order and index of 
g(mod p) over residue classes 

Pieter Moree 


Abstract 

For a fixed rational number g ^ {—1,0,1} and integers a and d we consider 
the set N g (a,d) of primes p for which the order of g(mod p) is congruent 
to a(mod d). For d = 4 and d = 3 we show that, under the Generalized 
Riemann Hypothesis (GRH), these sets have a natural density 6 g (a, d) and 
compute it. The results for d = 4 generalise earlier work by Chinen and 
Murata. The case d = 3 was apparently not considered before. 


1 Introduction 

Let g ^ { — 1,0,1} be a rational number (this assumption on g will be maintained 
throughout this paper). For u a rational number, let u p (u) denote the exponent 
of p in the canonical factorisation of u (throughout the letter p will be used 
to indicate prime numbers). If v p (g) = 0, then there exists a smallest positive 
integer k such that g k = l(modp). We put ord 9 (p) = k. This number is the 
(residual) order of g(modp). The index of the subgroup generated by g mod 
p inside the multiplicative group of residues mod p, [(Z /pZ) x : (^(modp))], is 
denoted by r g (p) and called the (residual) index mod p of g. Although ord s (p) 
and r g (p) satisfy the easy relation 

ord g (p)r g (p) =p- 1, (1) 

the functions themselves fluctuate quite irregularly. Given this it comes perhaps 
not as a surprise that a simple question such as Artin’s primitive root conjecture 
(1927), which asserts that {p : r g (p) = 1} is infinite if g is not a square, remains 
unsolved. On the assumption of the Riemann Hypothesis for a certain class of 
Dedekind zeta functions, however, this was proved by C. Hooley HU- Many vari¬ 
ations of Artin’s conjecture have been considered in the course of time, the most 
far reaching in H71 

Many authors studied the divisibility of the order by some prescribed integer 
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d. The case d — 2 for example is closely related to the non-divisiblity of cer¬ 
tain integer sequences by a prescribed prime. We say that an integer sequence 
S = {s : , }jl 1 is divisible by an integer m, if there exists an integer k such that m|sfc. 
It is easy to see that for a prime p with v p (g) = 0 the sequence S(g) = {gi + 1}°?^ 
is divisible by p if and only if ord 9 (p) is even. Hasse fH, TDj showed that the set 
of prime divisors of the aforementioned sequence has a Dirichlet density. It is 
not difficult to extend his argument to show that these sequences have a natural 
density (hereafter we merely write density instead of natural density) of prime 
divisors. For g — 10 a prime p 2,5 divides S(g) if and only if the period of the 
decimal expansion of 1/p is even, cf. |28| . Using some algebraic number theory 
these results can be extended to some other well-known sequences, cf. mmm- 
In all of these cases the density can be computed unconditionally and turns out 
to be a rational number. For example, the density of prime divisors of S( 2) is 
17/24. 

Now let d > 2 be given. By similar methods the divisibility of the order by 
d or the coprimality of the order with d can be studied. In this direction we 
especially like to mention K. Wiertelak, who wrote many papers on this subject, 
starting in the seventies of the previous century. See m for his most recent 
paper. Again one can prove that the density of the set of such primes exists and 
is rational. 

In the light of the extensive literature on the case where the order is divisible 
by d, it is somewhat surprising that the question of how the order is distributed 
over the various residue classes mod d has up to this century only been considered 
for d — 2. The purpose of this paper and its sequel(s) is to address this question 
for various other values of d. For the understanding of the general case it is in 
my viewpoint crucial to first study a particular case in detail, for which we take 
d = 4. 

For d — 4 our main interest is in the set N g (a, 4), but it turns out to 
be fruitful to consider N g (l, 2 s ; j, 4)(x) and N g (3,4; j,4)(x) separately, where 
N g (ai, dp, 02 , d 2 )(a;) counts the number of primes p < x satisfying v p (g) = 0 
for which p = ai(mod d\) and ord 9 (p) = a 2 (mod d 2 ). For convenience we de¬ 
note N g (0,l-,a,d)(x) by N g (a,d)(x). Although the functions N g (l, 2 s ; 1,4)(x) 
and N g (l, 2 s ; 3,4)(a;) are more complicated (see Theorem EJ) to describe, they 
turn out to be asymptotically equal under GRH. For the more easily describable 
functions N g ( 3,4;l,4)(x) and IV 9 (3,4; 3,4)(x) (vide Lemma fTTTl) . the asymptotic 
behaviour can be different. 

For s|r the number field Q(Cr,5' 1 ^) will be denoted by K rs . By til(x) we 
denote the number of rational primes p < x that are unramified in the number 
field L and split completely in L. As usual we let Li(x) denote the logarithmic 
integral, that is Li(x )=fj 
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For the primes p = l(mod 2 s ), s > 2, we find: 


Theorem 1 Write g = gi/g 2 with g\, r/ 2 integers. Let s > 2. For j = 0 and 
j = 2 we have 


N g ( 1, 2 s ; j, 4) (x) = 5 g (l, 2 s ; j, 4)Li(x) + O 


a;(logloga;) 
log 3 a; 


( 2 ) 


where 


S g ( 1,2 s ; 0,4) 




and 


( 3 ) 


S (l 2 s - 2 4) — I_I_1. 

9 r>s ' \^2 r , 2 r ~ 1 '■ Q] [-f^2 r +l,2 r - 1 : Q] [-K'2 r ,2 r : Q] [i^r+l^r : Q] ' 

For j = 1 and j = 3 we have, under GRH, 

where 

6 g (l,2 s -,l,2) = g { [1<2r 2r :Q] - [av +1)2 , Tqj) 

is the density of the set N g { 1, 2 s ; 1, 2) and the implied constant is absolute. 

For the primes p = 3 (mod 4) we find: 

Theorem 2 Write g — g\/g 2 with g\ and g- 2 integers. Let if o, if i denote the prin¬ 
cipal, respectively non-principal character mod 4. Let h,^ (v) = p{d)if\{v / d), 

where p denotes the Mobius function. We have 

f N g (3, 4; 0, 4)(x) = 0; 

{ N g (3,4;2,4)(x) = #{p <x:p = 3(mod 4), (J) = -1}. 

Assuming GRH we have, when j is odd, 

N g {3,4-j,4)(x) = <x :p = 3(mod 4), (-) = 1} 
z p 


where 


+ (-l) 


7-i x 


4 log a; 


+ 0 


A » = E 


ft*, (v) 


'/-2£ k 2 V ,2v 

2\v 


[R2v,2v ■ 


log 1.9 it/21 



5 


- E 

^^ K 2v,2v 

2\v 


[-^2d,2j; : Q. 


and the implied constant is absolute. 
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It is clear that N g (j,4)(x) = N g (l,A;j,A)(x) + N g (3, 4; j, 4)(x) and we leave it to 
the reader to add the estimates for the latter two quantities given in Theorem |T| 
respectively Theorem E] (Corollary U provides an example). 

In Section [21 we derive explicit versions of Theorem [7] and Theorem [21 For 
s > 0 it follows that 5 ff (l, 2 s ; j, 4) exists and is in Q + QA^, where 

A- = TT ( 1 - - % - 1 = 0.643650679662525 ... 

h 11 V p 2 + l P-1 / 

p=3(mod 4) V ’ 7 

As an example we mention the following corollary to Theorem [T| and Theorem [21 
(for the notation h and D we refer to Lemma El for a non-zero real number r we 
denote its sign by sgn(r)). 


Corollary 1 (GRH). Suppose that h = 1 and j is odd. Then h g (l,4) = d g (3,4) = 
1/6 unless D is divisible by 8 and has no prime divisor congruent to l(mod 4), 
in which case we have 


SgidA) 


2+1 

I + sgn(g)A^ Yl p] f jP-Xp- 1 tf D ± 8; 

3+1 

is + sgn (g) Afa if D = 8. 


In case d — 3 similar results to those for d = 4 are obtained in Section El In 
particular we will show that, under GRH, 5 9 (1, 3 s ; j, 3) exists for s > 0 and that 
Sg( 1, 3 s ;j, 3) G Q + QA a , where 


+. - n 

p=2(mod 3) 


2 P \ 

{. V 2 + i)(p- 1)/ 


0.173977122429634..., 


and denotes the non-principal character mod 3. The rational numbers involved 
we explicitly compute. 

The analogous problem of studying the primes for which the index is congru¬ 
ent to a(mod d) turns out to be far easier (see Section El). Nevertheless, at least 
for d = 3 and d = 4, we again find, under GRH, that these densities exist and 
are in Q + , respectively Q + . 

Instead of requiring GRH it is enough to require that RH holds for every held 
Q(Cr,# 1 ^) with s\r. Indeed, if a given result is under GRH we mean that we 
require RH to hold for every held that occurs in the proof of this result. 

In a sequel to this paper by a slightly different method the case where d is an 
odd prime power is investigated (but less explicitly), see |23j. 

The density 5 g (j, 4) with g a positive integer that is not a pure power (i.e. 
h — 1 in the notation of Lemma [TJ , was hrst studied by Chinen and Murata in 
El II El 0 0 EZI, culminating (in EH) in their proof of Corollary 1 for the case 

g> o. 


2 Preliminaries 

2.1 The index and algebraic number theory 

In this section we recall some well-known arguments from the theory of primitive 
roots that are essential for an understanding of the rest of this paper. 
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The index can be easily related to algebraic number theory and by using © 
we then can get a grip on the order. Thus a prime p that satisfies k\r g (p) must 
obviously satisfy p = l(mod k) and gG~ 1 )/ k = l(mod p), in other words it must 
split completely in the field Q{Ck, g 1 ^ k )(= Kk,k)- On the other hand a prime 
p that satisfies the latter condition satisfies k\r g (p). Then, by the principle of 
inclusion and exclusion, we can describe for example the set of primes p that 
satisfy r g (p ) = k. Note that r g (p ) = k iff k\r g {p) and qk \ r g (p ) for any prime 
q. Let R g (a, /; t ) denote the set of primes p with p = a(mod /) and r g (p ) = t. 
Let R g (a, f]t)(x) denote the number of primes p < x in R g (a, f',t). Using the 
principle of inclusion and exclusion, we then find that 

OO 

Rg(a,f]t)(x) = ^/fyn)#{p < x : p = a(mod /), (p, K ntM /Q) = id}. (4) 

n= 1 

By (p,K/ Q) we denote the Frobenius symbol. We have (p,K/ Q) = id iff p 
is unramified and splits completely in K. Since sets of the form {p : p = 
a(mod /), (p, K rtn /Q) = id}, will occur rather frequently in the sequel, we 
will denote them by S g (a, f',r,n) and the corresponding counting function by 
S g (a, f;r,n)(x). Assuming GRH, it follows from ;!?] that R g (a, f;t ) has a den¬ 
sity. 

Sofar this density has only been evaluated in terms of an Euler product (sin¬ 
gular series) in the case t — 1 with a and / arbitrary [TSj, or in the case f\2 and 
t arbitrary mm- For example, for t = 1 and 2\f the density, under GRH, is a 
rational multiple of the Artin constant f([ p (l — ) [11 j . For an unified Galois 

theoretic treatment of finding Euler products for these cases see m- 

For our purposes such an evaluation of the density of R g (a, /; t) will, however, 
be irrelevant, an evaluation in terms of an infinite series will be sufficient. The 
tool to arrive at such an expression for the density is the Chebotarev density 
theorem: 

Theorem 3 (GRH). Let K be an algebraic number field, let L/K be a finite Ga¬ 
lois extension and C be a conjugacy class inG = Gal (L/K). We let 7rfy; L/Ii, C) 
denote the number of unramified prime ideals p in K such that (p, L/K ) = C and 
Np < x. Then, wider RH for the field L we have 

tt(x-,L/K,C) = J^yLify) +0 ^^V^log (d L x [L:Q] )^j , as x -> oo, (5) 
where dL denotes the discriminant of L. 

Remark. The proof of Theorem [3 is in essence due to Lagarias and Odlyzko, the 
present formulation is due to Serre m p. 133], who removed ‘un terme parasite’ 
in the formulation of Lagarias and Odlyzko (HSl Theorem 1.1]). In case C = id 
the result was proved earlier by Lang M- For several variants of Artin’s primi¬ 
tive root conjecture Lang’s result is all one needs. There are also unconditional 
variants that certainly allow us to deduce that n(x; L/K, C) ~ 4§Li(x), as x 
tends to infinity. 
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For an arbitrary integer m > 1 let us see how Chebotarev’s density theorem 
can be used to estimate S g (a, f;r,n)(x), where n\r. To this end we consider 
the compositum of the fields Q(C/) and K r ,n, that is where by [/, r] we 

denote the lowest common multiple of / and r. Let K\ and K 2 be number 
fields that are Galois. If there is an automorphism oq G Gal(A'i/Q) and an 
automorphism <j 2 £ Gal(iL 2 /Q) such that = o 2 on K\ fl K 2 , then there is 
an unique a G K 1 ■ K 2 , the compositum of K\ and K 2 such that a\x = &j for 
j = 1,2. Now in order to apply the Chebotarev density theorem, we have to count 
the number of elements in the conjugacy class of a G Gal(AT[/ )r .] jn /Q), where a 
is such that ct|q(q) = cr a j, where <r a j G Gal(Q(C/)/Q) is uniquely determined 
by (Taj((f) = (f, and a a j \ j< rn = id. By the above remark such a a exists, and 
is unique, if and only if cr|Q(C/)n K r , n = id- Note that a conjugate rar -1 acts 
trivially on K r n and can be regarded as an element of Gdl(Kj :r i n /K r)n ), which is 
a subgroup of the abelian group (Z//Z)*. Hence rar -1 acts as a. We conclude 
that the conjugacy class has one element if 0'|Q(c / )niq.,„ = id and zero otherwise. 
By this argument we expect from 0, assuming there is enough cancellation in 
the error terms, that the density of R g (a, /; t ) is given by 



where 



and n(x) = J^ P < X 1 - % HU we know that © holds true, under GRH. 

2.2 Field degrees and intersections 

In order to explicitly evaluate certain densities in this paper, the following result 
will play a crucial role. The notations D, g 0 , h and n r will reappear again and 
again in the sequel. If a and b are integers, then by (a, b ) and [a, b] we denote the 
greatest common divisor, respectively lowest common multiple of a and b. 

Lemma 1 Write g = ±<7o., where g 0 is positive and not an exact power of a 
rational. Let D denote the discriminant of the field Q( v /^o). Put m = D/2 


if v 2 (h) = 0 and D = 4(mod 8) or u 2 (h) = 1 and D = 0(mod 8), and m = 
[2 U2 ^ h ) +2 , D] otherwise. Put 


m 


if g < 0 and r is odd; 


[ 2 V 2 (.hr)+i ^ d] otherwise. 


We have 



where, for g > 0 or g < 0 and r even we have 
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and for g < 0 and r odd we have 

( 2 if n r \kr; 

e(kr, k) = < 1 if 21 k and 2 U2 ^ +1 \ k; 

{ 1 otherwise. 

Proof For r = 1 the result follows from Proposition 4.1 of jHJ (see also the proof 
of Theorem 2.2 in EH). For r > 1 the result follows from the case where r = 1 
on noting that, with g = g r , we have Q(Cfcr, g 1//fc ) = Q(Cfcr, g 1 ^ kr ). The distinction 
between r even and r odd arises in the case g < 0 since g is then positive or 
negative, according to whether r is even or odd, respectively. □ 

In our analytic considerations we need an upper bound for the discriminant 
of the held K krk . 

Lemma 2 The discriminant D' of the field K r k : k satisfies 

log \D'\ < rk(\og{rk) + log(|^ 2 |)), 
where g 0 = g k /g 2 and g\ and g 2 are integers. 

Proof. If L i/Q and L 2 /Q are two extension fields and L is their compositum, 
then the associated discriminants (over Q) satisfy d[fd^ L '\j\ L 2 L2 K From this we 
have the estimate 


log Ml I < [L 2 ■ Q] logical + [Li : Q]log|d i2 |. (7) 

It is well-known that the discriminant of the cyclotomic held Q(Cm) and the held 
Q(fi ,1 ^ n ) divide m respectively ( ngig 2 ) n (see e.g. j21J). On invoking these es¬ 
timates, the result then follows from O □ 

From cyclotomy we recall the following well-known result. 

Lemma 3 We have Q(y/g) C Q(C/) iff A|/ ; where A denotes the discriminant 
of the field Q (y/g). 

Proof. See e.g. |J2]. □ 


In order to evaluate the densities for the modulus 4 we need the following 
result, which can be easily deduced from Lemma H] and the previous lemma. 


Lemma 4 Put L v = Q(Cs) 0 Q(C 2 v, g 1 ^ 21 ’)■ Let v be odd. 
If h is odd and D \ 8v, then L v = Q. 

If h is odd and _D|8u, then 


If h is even, then L 


Q(\ /sgn (g) ) 
Q(\ /sgn (-g) ) 
Q(v / 2-sgn(^) ) 
Q(v / 2 T sgn (Pg) ) 

= Q(y / sgn (g) ). 


if D = l(mod 4); 
if D = 4(mod 8); 
if D = 8(mod 32); 
if D = 24(mod 32) 
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Proof. On using that 


[Q(Cs.,^ 1/2 ") 


4 


mC2v,g 1/2v ): Q] 

[L v : Q] 




it follows by Lemma[T]that [L v : Q] = e(2v, A)/e(2v, 1). Using Lemma Q] it follows 
that in each case the field claimed to equal L v has the correct degree. 

Let us first consider the case where g > 0. We may suppose that h is odd 
and D\8v, since in the remaining cases the degree of L v is 1 and hence L v = Q. 
It remains to show that i G L v if D = 4(mod 8), \[2 G L v if D = 8(mod 32) and 
v 7 —2 G L v if D = 24(mod 32). 

i) D = 4(mod 8). In this case the discriminant of Q(y/—g), D/A , divides v and 
thus by Lemma 0 y/—g G Q(Cu) Q K 2Vj2v . Since also y/ g G K 2Vj 2 v it follows that 
i G K 2v 2v . 

ii) D = 8(mod 32). Now a/2 ~g G Q(( v ) and thus \/2 G K- 2V)2v . 

iii) D = 24(mod 32). Now a/— 2 g G Q(C„) and thus a/^2 G K 2v2v . 

Suppose g < 0. If h is even, we have to show that i G Q(C 2 «, g 1 ^^) = 
Q(C 2 vX 4 v% / 2v )- Since (Cav9q^ 2v ) v G Q(*)\Q, this is clear. As before we may now 
suppose that h is odd and D\8v. It remains to show that i G L v if D = l(mod 4), 
V—2 G L v if D = 8(mod 32) and y/2 G L v if D = 24(mod 32). Note that 
(( 4 v% /2v ) v is a rational multiple of y/— g 0 and hence a /—go G K 2v2v . 

i) D = l(mod 4). We have y/go G Q(Cu) and since y/—go G K 2Vj2v , it follows that 
i G K 2v 2v 

ii) D = 8(mod 32). We have y/2go G Q(Cu) and since a/— ~go G K 2v , 2v , it follows 
that a/—2 G L v . 

iii) D = 24(mod 32). We have \/—2g 0 G Q(</) and since a/— G K 2V)2in it 

follows that y/2 Gl„. □ 


Lemma H] allows one to establish the following property of A g . 

Proposition 1 We have A_ g = —A g . 

Proof. If h is even, then A_ 9 = A 9 = 0. If h and v are odd and 8|.D then 
[Q(C 2 ^ 1/2 " : Q] = [Q(C 2 u, (— g) 1 ^ ■ Q] by Lemma [U The result now follows 
easily on invoking Lemma |U □ 


2.3 Index t revisited 


In this section we extend some results of Murata ESI, which he established for 
squarefree integers > 2, to arbitrary g G Q\{ —1,0,1}. However, the method of 
proof employed in this paper is rather different. 

We will need that ^ ~f r>x l/(rip(r)) = 0(1/x) (for a proof see e.g. O Lemma 
8.4]), the following result shows that this order of growth cannot be improved 
upon. 


Lemma 5 We have 


£ 

n<x 


l 

mp(n ) 



P \ 
1 


315C(3) /logs 

2n A x V x 2 














Proof. Landau na has shown that 


E 

n<x 


1 

<p{n) 


315C(3) 

2vr 4 


{logx + 7 - 


E 


logp 

p 2 — p + 1 


} + 0 



( 8 ) 


where 7 denotes Euler’s constant. Using © and the Euler identity, the result 
follows on partial integration. Alternatively one can apply Landau’s method for 
establishing (JHJ) to the sum J2n<x 1 /(n<^(n)). □ 


The estimate 


ensures that 



< 


2 h 

vtip(vt ) ’ 


A(g,t) 


E 

V=1 


[Kvt,vt '■ Q] 


(9) 


converges absolutely. 

Lemma 6 We have Y2t<y A(g,t) = 1 + O(^), where the implied constant is 
absolute. 


Proof. One easily checks that Y!itLi-A-{g,t) = L Using 0, <f(vt) > ip{v)ip{t) and 
Lemma El we infer that 


EE 


h 


t>y 


< v 


Vt<p(vt) 


o (v—— v — 

1 Mt) ^ v<p(v) t 


= o(~), 
y 


where the implied constants are all absolute. 


□ 


Theorem 4 (GRH). Write g = gi/g 2 , with g± and g 2 integers. Fort < x 1 / 3 
have 


R g (0, l;t)(z) = A(g,t)—^— + 0 

logx 

where the implied constant is absolute. 


f X log logx glog|glff 2 | \ 
V <p(t) log 2 X log 2 X ) 


1 


we 


Proof. Since the proof is carried out along the lines of Hooley’s proof m , we 
only sketch it. Let 

M g (x,y) = #{p < x : t\r g (p), qt\r g (p ), q < y} and 

M g (x, y , z) = #{p < x : qt\r g (p), y <q< z}, 
where q denotes a prime number. Note that 

R g ( 0,1 ;t)(x) = Mg(x,T 1 )+0(M g (x,T 1 ,T 2 ))+0(Mg(x,T 2 ,T 3 ))-l-0(M g (x,T 3 , —j—))- 

We take 17 = logx/ 6 , r 2 = \/xlog~ 2 x and r 3 = yfx logx. We use the starting 
observation that 

OO 

R g (0,l]t)(x) = ^2y{v)n KvtM (x). 

V=1 
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Using Lemma [2] and Theorem 01 we can estimate 7 TK vt , v t( x )i under GRH. Proceed¬ 
ing as Hooley did, we then obtain that both M g (x, £ 1 , C 2 ) and M g (x, are of 

order (log \gig 2 \)x log -2 x, where the implied constant is absolute. Furthermore, 
we obtain that M g (x, ( 2 , C 3 ) — w h ere again the implied constant is 

absolute. For the main term we find that 


M g {x, Ci) = A(g,t)-^— + O 
logx 


{ ^ log |gl^2l A 
V log 2 X ) 


where the implied constant is absolute. On adding the various terms, the theorem 
follows. □ 


The following result is a slight generalisation of Lemma 2.4 of ITj. 

Lemma 7 (GRH). Let 'ijj(x) be a monotonous increasing positive function which 
satisfies 

lim = +00 and ip(x) -C (logx) 1 / 2 . 

X—>00 

Then we have 

nix) 

#{.P < X : r g {p) > if{x)} < log \9i92\^tt, 

where the constants implied by the <C -symbol are absolute and g — gi/g-i with g\ 
and g 2 integers. 

Proof Let y denote the largest integer not exceeding -0(x). We have 

#{p < x : r g (p) > y} = {p < x : v p {g) = 0} - U^'P^O, 1 ;t)(x), 

where U^L 1 is a disjoint union. The latter identity together with Theorem 0] 
Lemma ED and (0), yields 

#{p < x : r g (p) >y} = tt(x) + 0(log \gig 2 \) 

K\ X + 0 f x logy log log x \ + xy log I 1 \ 
y J logx V log 2 x J V log 2 x J ’ 

where the implied constants are all absolute. On using the prime number theorem 
in the form 7r(x) = xj logx + 0(x log~ 3 / 2 x), the result then follows. □ 



2.4 On the convolution of the Mobius function with Dirich- 
let characters 

Let y be a Dirichlet character of conductor f x and order o x (for definitions and 
basic facts on Dirichlet characters we refer the reader to Hasse 0). Let Gd denote 
the group of characters defined on (Z/oK)*. We have that Gd — (Z/oK)*. An 
important auxiliary function in this paper is the convolution, h x = y* x, of the 
Mobius function p with a Dirichlet character x, that is h x (n) = J2d\n p(d)x( n /d). 
In this section we collect some auxiliary results involving h x . We note the fol¬ 
lowing trivial result. 
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Lemma 8 The function h x is multiplicative. With the convention that 0° = 1, 
it satisfies h x { 1) = 1 and h x (p r ) = x{p) r ~ l {x{p) ~ !]• 

In particular if \ is the trivial character mod d, then 

h (v) = I if v\d\ 

lO otherwise. 

By using one of the orthogonality relations for Dirichlet characters, the following 
result is easily obtained. 

Lemma 9 Let a(mod d) be a reduced residue class mod d. We have 

Y, = X(a)h x (v), 

t=o(mod d) T\ ) x (iQ d 

t\v 

where Xk runs over the Dirichlet characters modulo d. 


Note that the lemma expresses a non-multiplicative function as a linear combi¬ 
nation of multiplicative functions. This will play an important role later on. 

Let r, s be non-negative integers. Put 


C x (h, r, s) 


E 


(r,t;)=l, s\v 


h x (v)(h,v) 

V(p(v) 


and A x — 


x(p)^0 



[x(p) ~Mp \ 
\p 2 - xip)\ip- i)y ’ 


It is easy to see that the latter series is absolutely convergent. Note that h x (v) < 
2 U ( V \ where u(y) denotes the number of distinct prime divisors of v. Note also 
that for every e > 0 we have \h x (v)\ < < Yld\v 1 <C v e and (p(v) 3> u 1_e . 

From this the absolute convergence easily follows. Since h x is a multiplicative 
function, we can invoke Euler’s identity. After some tedious but easy calculations 
this then yields the following result. 


Lemma 10 Let h, r,s > 1 be integers. Let x be a Dirichlet character mod d. 
Put e p = v p {h). 
i) If (r, s) — 1, then 


C x (h, r, s) = 

p\rs 


l + p 1 ep 


\x(p) - 1] 
ip- 1) 


f p ep - x{p ep ) x(p ep ) \ \ 
V p-xip) p 2 -x(p)JJ 


n p 1 ep 

!/p(s) = l 


[x(p) - !] 
ip~ !) 


/ p e P _ x ( p e p ^ x(p £ p) \ 

V p-xip) p 2 -x{p)J 


n 

vp(s)>ep-\-i 

vp(s)>2 




[x(.P) — 1] p e P+3 2vp(s) 

ip- 1 ) (p 2 ~x(p)) 


tt [x(p) — i] f {x{p)/pY 2[s) l -{x{p)/p) ep | p 1 ep x{p ep ) \ 
"d W ip - x ) v 1 - xip)/p P 2 - xip) ) ’ 


In particular C x (h, r, s ) = cA x = cC x ( 1, d, 1), with c G Q(Co x )- 
ii) If (r, s) > 1, then C x (h, r, s) — 0. 
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Remark. In this paper we only need to evaluate C x (h,r,s) in the case where the 
largest odd divisor of s is squarefree and ^(s) > e 2 + 1, in which case it assumes 
a bit simpler form than the general one given in Lemma Hill 


Only the primes p with x(p) 7^ 1 contribute to A x . Note that A x ^ 0. If x is 
a principal character, then A x = 1. If x is real, then A x > 0. If x is a Dirichlet 
character and x! is its associated primitive Dirichlet character, then clearly A x 
equals A x i with some local factors involving only ( 0x left out. Hence A x = c'A x > 
with d G Q(Co x )- Also note that A x = A x . 

The constants A x are the basic constants in this paper. They have a product 
expansion in terms of special values of Dirichlet L-series j2HG2J. This result is 
related to the denominator identities arising in the theory of Lie algebras 1231 . 
These expansions can be used to evaluate A x with high numerical accuracy I2Tj . 
The values of A^ and A^ 1 given in the introduction are taken from Table 3 of 

Ea- 

Another result involving h x needed is the following (where the sums are over 
the integers v > 1). 


Lemma 11 Let r, s be integers with s \r and g > 0. Let x be a Dirichlet character. 
Then, if g > 0 or g < 0 and s is even, 


h x (v) 

[K sv v : Q] 

(r,v )=1 


<p(s) 


C x (h,r ,!) + C x (h,r, 


n. 


[ n si s ) ’ 


When g < 0 and s is odd, the latter sum equals 


1 

<P(s) 


C x (h,r,l ) 


tc x (h,r,2) + l -C x (h,r,2^+ 1 )+C x (h,r. 


-2i_)V 

K,s) ) 


Proof. For g > 0 or g < 0 and 2|s, the proof easily follows from the identity 


E 

(r,u)=l 



E 

(r,«)=l 


h x (v)(h,v) 

vtp(sv) 


+ E 

(r,v) = 1 
n s \sv 


hx(v)(h,v) 

inp(sv) 


which on its turn is an easy consequence of Lemma |U The proof of the remaining 
case is similar. □ 


2.5 Preliminaries specific to the case d = 4 

2.5.1 Even order 

Let s > 2. It is not difficult to estimate IV 9 (1,2 s ; 0, A)(x). To this end we 
consider the set of primes p = l(mod 2 s ) such that ord 3 (p) ^ 0(mod 4). Let 
r = V 2 (p — 1). Note that r > s. Now p satisfies ord 5 (p) ^ 0(mod 4) if and only 
if g ( -P~ l T 2 ' r = l(modp), that is if and only if p splits completely in AVg 1 - 1 , 
but not completely in K 2 r+i^ r ~ 1 - By Chebotarev’s density theorem we expect 
that IV 9 (1,2 s ; 0,4) has a density as given in (|3|). Indeed, this can be shown 
unconditionally. Wiertelak’s work [ 331 Theorem 2] goes beyond this and shows 
that we even have an estimate for iV 9 (l, 2 s ; 0,4)(a;) as given by P|). 

Similarly we can easily estimate N g {l, 2 s ; 1,4)(a;) + iV 9 (l, 2 s ; 3,4)(x) which 
equals IV 9 (1,2 S ; l,2)(x) and AT 9 (1, 2 s ; 2,4)(x). See Theorem[I] for the outcome. 
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2.5.2 Odd order 


We let ipo, denote the trivial respectively non-trivial character mod 4. The 
starting point of our analysis is the following easy result. 

Lemma 12 Let s > 1. For j odd we have 
N g (l,2 s ]j,4:)(,x) = ^2 < X: P = 1 + 2 r (mod 2 r+2 ), r g (p) = 2 r t} 

r>s t=j (mod 4) 

+ ^2 ^2 < x : p = 1 + 3 • 2 r (mod 2 r+2 ), r g (p) = 2 r t}, 

r>s t=—j (mod 4) 

N g (3,8;j,A)(x) = ^ #{p < x : p = 3(mod 8), r g (p) = 2 1} and 

t=j (mod 4) 

Ng(7,8]j,4:)(x)= ~22 #{p < x :p = 7(mod 8), r g {p) = 2t}. 

t=—j (mod 4) 

Proof. We only prove the assertion regarding iV s (l, 2 s ; j, 4)(x), the other asser¬ 
tions being easier to prove. For every prime p = l(mod 2 s ) there exists an unique 
r > s such that either p = 1 + 2 r (mod 2 r+2 ) or p = 1 + 3 • 2 r (mod 2 r+2 ). We 
assume that we are in the first case, the other case being dealt with similarly. 
Using 0 we note that ord g (p) = j(mod 4) if and only if r g (p) = 2 r t for some 
t> 1 with t = j (mod 4). □ 

The densities of the sets appearing in Lemma [12] can be determined, under GRH, 
on invoking Assuming the densities add up, we then arrive at the conjecture 
that the densities are as stated in Theorem |H1 (cf. Section HJ). 

An alternative approach is obtained on first resumming the expressions in 
Lemma d and only then applying Lenstra’s machinery mi. which is what is 
explained next. 

Let s > 1. From Lemma d and m we deduce that, when j is odd, 

OO 

N g (l,2’-,j,4)(x)= Y. yjMn)S 9 (l + 2',2’'+ 2 in«2',nt2')W 

r>s t=j (mod 4) n= 1 


oo 

+ S 5Z J2p{n)S g {l + 3-2 r ,2 r+2 ;nt2 r ,nt2 r )(x). (10) 

r>s t=—j {mod 4) n= 1 

Note that the fields that arise in the S g occurring in cdd are of the form K V 2 r , V 2 r - 
On grouping together the contributions involving the various K v2 r, V 2 r the triple 
sums can be reduced to double sums. To this end we first note that we can 
restrict to the case where n is odd, since if n is even and m is odd, S^l + 
m ■ 2 r , 2 r+2 ; nt2 r , nt2 r ) is empty (then the condition on the Frobenius symbol 
implies that p = l(mod 2 r+1 )). On putting v = nt the summation is then over 
all odd v > 1. As weighing factors we then get sums as in Lemma El with 
(a,/) = (1,4) and (a,/) = (3,4). On applying Lemma El and noting that for 
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odd v, h^ 0 (v) — 1 if v — 1 and h^ 0 (v) = 0 otherwise, we then obtain that 
N g (l,2 s ]j,A)(x) = |/i + ( ^ 1} 2 ^~ / 2 , where 

h = ^ [#{p < £ : (p, AV,2-/Q) = - #{P < a; : (p, K 2 r+ i i2 r/Q) = id} 

r>s 


and 


= V ( v ) [^(! + 2 r > 2 r+2 ; y 2 r , ^2 r ) (x) - (1 + 3 • 2 r , 2 r+2 ; AT, vT) (x) 

r>s 2\v 


(Note that the latter double sum can be simplified to a single sum. On doing so 
we find that I 2 equals 


£ 


S s (l+2 ra (">,2«<»> +2 ; ro , ro )(i)-S ! ,(l+3-2 ra <»>,2"(”)+ 2 ; ro , ro )(i) 


where w 0 dd is the largest odd divisor of w.) Note that if we add N g (l, 2 s ; 1,4)(x) 
and N g (l, 2 s ; 3,4)(x) we obtain J 1; which is a well-known result. 

If s > 2, then the Chebotarev density theorem implies, unconditionally, that 
for v odd and r > s, 

S g { 1 + 2 r , 2 r+2 ; u2 r , u2 r ) (x) ~ S g ( 1 + 3 • 2 r , 2 r+2 ; v2 r , v2 r ) (x), as x -»• cx). (11) 


Thus we might expect that I 2 behaves like an error term and that, consequently, 
7V s (l,2 s ;l,4)(x) ~ I\/2 as x tends to infinity. Theorem ^ shows that this is 
indeed true, under GRH. 

If s = 1, however, then m does not necessarily hold true for every r > s. 
It thus makes sense to consider N g (3, 4; j, 4)(x) for j odd separately. We then 
obtain analogous expressions to those for I\ and I 2 , but instead of summing over 
r > s we take r to equal one: 

N g (3, 4; j, 4)(x) = ^#{p <x:p = 3(mod 4), (-) = 1}+ 
z p 


i\ti 
- 1 ) 2 


v ( v ) 

2\v 


S g (3, 8; 2v, 2v)(x) — S g (7, 8; 2v, 2v)(x) 


Using that p splits completely in Q(v / —2) iff p = l(mod 8) or p = 3(mod 8) and 
that p splits completely in Q(y / 2) iff p = ±l(mod 8), we obtain the following 
result. We present it with a more succinct proof. For a number field L we let 
7Tt,(x) denote the number of rational primes p < x that split completely in L. 


Lemma 13 Let j be odd. For every x we have 


As(3,4;j,4)(x) = < x : p = 3(mod 4), (^) = l}d 2 


1 

1) 2 


^h in {v)A v (x), 


2\v 


where A v {x) = tt K2V}2v (V^)( x ) ~ *k 2v mV 2 )( x )- 
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Proof. Let us consider only the case where j = l(mod 4), the remaining case being 
dealt with similarly. Note that in both the left hand side and the right hand side 
of the identity that is to be established only primes p satisfying p = 3 (mod 4) and 
v p (g) = 0 are counted. Now let p be a prime such that p = 3(mod 4), v p (g) = 0 
and p < x. We shall show that it is counted with the same multiplicity in both 
the left and the right hand side of the identity that is to be established, thus 
finishing the proof. 

Since by assumption v p (g) = 0, there exists a largest integer k such that 
g~k~ = l(mod p). Note that in both the left and the right hand side only primes 
p with k even are counted. Thus we may write k = 2 k\. Note that k\ must be 
odd. Let us assume that p = 3(mod 8). Then p is counted on the right hand side 
with weight 


1 1 V - -' , ,N 1 l /( i\/r\ l + ^l(fcl) 

2 + 2LVH = g + 2* R * x )(^i) =-2-' 

m\ki 

Thus the weight is 1 if k\ = l(mod 4) and 0 otherwise. In other words the weight 
is 1 iff ord 9 (p) = l(mod 4). 

The case where p = 7(mod 8) is dealt with similarly. □ 


Using Chebotarev’s density theorem we now expect that 

(u) hpi (u) 


Hm h^(v)A v {x) _ ^ 


n(x) 


2\v 


[A 2v,2v(\ /—2) : Q] [K2v,2v(V 2) : 


Let us denote the quantity on the right hand side by A„. Note that 


A, = - 
9 2 


E 


V zz 2eK 2v ,2v 

2\v 


h 1 (t>) 

[k^2v,2v '■ ^ 


E 


K Av) 


A, 


2 [P2v,2v '■ 

v2€ k 2v,2v 
2\v 


Theorem [21 shows that this heuristic holds true, under GRH. 


3 The distribution of the index over residue classes 


The problem of the distribution of the index over residue classes is far easier than 
that of the distribution of the order. However, the answers to both problems turn 
out to have some features in common. 

Let a and d be integers. Under GRH it follows from Pappalardi’s work 
that the density, p g (a, d ), of the set of primes p such that r g (p ) = a(mod d) exists 
and equals 


Pg ((X, d) 


E E 

t=a(mod d) v=l 


p(v) 

[■ Kvt,vt ■ Q] 


Using this the following result is then easily deduced. 
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Theorem 5 (GRH). Let a and d be arbitrary natural numbers. Put 5 = d/(a, d). 
Then the density of the primes p with r g (p) = a(mod d), p g (a,d), exists and 
satisfies 

p g (a,d) = Y c x A x, with c x e Q(CoJ- 

X^Gg 

Furthermore, c x = c x . The number c x can be explicitly computed. 

Proof. On putting w = (a, d) and a = a/w we obtain 


Pg( a i d) — Y zL 


p(v) 


L vwt.vwt • 


, , £ \ . \Ky 

t=a( mod 5) v=l 

Writing vt = V\ and invoking Lemma 0 we then obtain 


p g (a,d) = Y X{<*) Y 


h x (v i) 


tp(6) 

Let g > 0. By Lemma [T| we have 


X&G S v-i 




VlW,ViW • 


E 


h x (v) 


L vw.vw • 


ti 1 A '« 

We can rewrite J\ as 


E 


h x (v)(h, vw ) 


E 


h x (v)(h, vw) 


—' VW(f(vw) vwcp(vw ) 

n\ | vw 


— J\ + J 2 


Ji = 


( h,w ) h x (v)(h, vw)<p(w) 


E 


wp(w)‘Y (h,w)v<p(vw) 


where the argument of the sum is easily seen to be multiplicative in v. If p \ hwf x , 
then the local factor at p in the Euler product for J\ equals that of A x and so 
Ji = c' x A x , where c' x £ Q (Co x )- Rewriting the condition nf\vw as |n, we see 
that also J 2 = c' x A x , where c" £ Q(Co x )- A similar argument can be used in the 
case where g < 0, cf. the proof of Le mm a fill 

On using that h x = h x and A x = A x it follows that = c x . □ 


A special case occurs when d\a. Then p g (a,d ) is the density of primes p < x 
such that the index of (/(mod p) is divisible by d, that is p g (a, d) is the density of 
primes p < x that split completely in K d d . By an unconditional version of the 
Chebotarev density theorem we then infer 

Proposition 2 The density p g (0,d ) exists and satisfies 

In the three examples some special cases of Theorem 0 are discussed. 

Example 1. (GRH). We consider the case where (a,d) = 1 and g > 0. Then, 
using Lemma ITTl we obtain that 

Pgifli d) 7 ,, E X(a) (c x (h , 1,1) + C x (h, 1, m)) . 

^ a) xe G d 
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Example 2. (GRH). We assume that ( a,d ) = 1 and g — 2. Note that D — 8 
and hence n± = 8. On invoking the formula of Example 1 with h — 1 and ri\ = 8, 
we obtain, using Lemma [TI71 


p 2 (a J cl) 


1 

ip(d) 


X&Gd 


with 


/I X<2)((x(2)) 2 - x(2) + 12)\ T-r / i>(x(p) - 1) \ 

8(4-x(2)) MU (P ~ 1)(P 2 ~ X(P))) ' 


(This corrects a typo in Corollary 8 of |29|). Alternatively we can write 


p 2 {a,d) 



P(P ~ !) 


*( a ) ( 1 + 

X'SGd 


X(4)(x(2)-1) \ 

8(x(2) + 2) J 


A x . 


Example 3 (GRH). Let g — 2 and d — 3. Using Example 2 we compute 
p 2 (±l, 3) = ^ ± that is p 2 (l, 3) = 0.471034 ■ • • and p 2 ( 2, 3) = 0.362298 ■ • ■. 
(This corrects the values claimed in p. 386].) From the latter two values 
we infer that (0, 3) = 1/6. By Proposition [21 it follows that we even have 
unconditionally that <5 2 (0, 3) = 1/6. Up to p 10 e = 1299709 the approximations 
0.16589, 0.47127 and 0.36283 for p 2 (0, 3), p 2 (l,3) respectively p 2 (2, 3) are found. 


4 Proof of the main results 

In our formulation of the main results stated in the introduction, we have added 
the (known) case d = 2 for completeness. As already pointed out in Section 
12 m these cases have been well studied. The best known error terms are due 
to Wiertelak, cf. [31] • The densities for these cases can be explicitly evaluated 
using Lemma El Since the interested reader can easily carry this out herself, we 
abstained from writing down the rather lengthy (because of case distinctions) 
outcome. For some further elaboration on these cases see Section 2.5.1. 

Our proof for the remaining cases has an analytic and algebraic component, 
with the analytic component being captured by the following result. 

Theorem 6 (GRH). Let s > 1. Let ij }i be the non-principal character modulo 
4. For r > 1 let cr_i )T . G Gal(Q(</ 2 r+ 2 )/Q) be the automorphisms that are 
uniquely determined by cr lir (C 2 H- 2 ) = C 2 ^+2 , respectively cr_ lir (C 2 r + 2 ) = 

For j = —1 and j = 1 let 

Cj (r, tn) = ( 1 ^ a P r Wc 2 r+ 2 )nQ(<: 2 rt n ,g 1/2rtn ) 

V 0 otherwise. 

For j = 1 and j = 3 we have 

N g {l,2 s -j,4)(x) = <5 g (l, 2 s ; j, 4)—^ f O (log\g 1 g 2 \-^- \ , 

logo; V log 6/2 xj 
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where 


«i.2';^4) = £££ 


H(n) c rhw ( r ’ tn ) 


r>s £= l 
— 2fi 


- [Q^C^n,^ 2 ^) 


and £/ie implied constant is absolute. 


An heuristic argument in favour of the truth of the latter theorem is easily given. 
From Lemma M and (J4|) we deduce that iV s (l, 2 s ; j, 4) (re) equals 

OO 

E E E At( w )#{p < x : p = 1 + 2 r (mod 2 r+2 ), (p, Kywtn/Q) = id}+ 

r>s t=j (mod 4) n=l 

OO 

E E E M n )#(P <x:p = l + 3-2 r (mod 2 r+2 ), (p, K 2 r t n,vtnlQ) = id}. 

r>s t=—j (mod 4) n =1 

The density of the inner sums is given, under GRH, by (ffij). Assuming the densi¬ 
ties add up and there is sufficiently cancellation in the error terms, we then arrive 
at the heuristic that the densities should be as claimed. Our proof of Theorem |U] 
is in the same spirit: 


Proof of Theorem |bl Let us denote the triple sum in the formulation of the result 
^ E E E' All constants implied by the O-symbols in this proof will be absolute. 
The first formula of Lemma can be more compactly written as 

Ng( 1, 2 s ', j, 4) Or) = E E + (2 - MjtW, 2 r+2 , 2 r t)(x). 

r>s 2ft 

On retaining only the primes with r g (p ) < y, we obtain 

N g ( 1,2 s ; j,4)(x) = T^y) + 0(#{p < x : r g (p) > y}), (12) 


where _ 

fiMAE R 1 (l + (2-Mjt))2 r X +1 Xt)(x). 

r>s 2 \t 

2 r t<y 

The function R g ( 1 + (2 — ifi(jt))2 r , 2 r+2 , 2 r t)(x) can be estimated in the same 
way as R g ( 0, l;£)(x) (see the proof of Theorem HJ). We find that, under GRH, 


R,( 1 + (2 - Mt))2 r , 2 r+2 ,2 r t)(x) = U(x) y(n)^ {jt) (r,tn) 

+0 (log(|0i0 2 |) 




n=1 i-^\s2 r + 2 , (2 r tn, g 1 R Ttri ) : 

X log log X 


log 2 x J + \<p(2 r t) log 2 x 


This, when substituted in m , yields 

™ (^) +0 ( Ilos , g g , logI ) + ° ■ 
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where we used (0 and 


y y y_-_ 

^ ^ 2 r trup(2 r tn) 

r>s 2 \t n =1 ' v 7 

2 r t>y 


h 


o yy- 

' mmx(mn) 

\m>y n =1 v 7 


oh 

y 


cf. the proof of Lemma El On taking y = \J\og x in Ha, the result then follows 
on invoking Lemma 0 with f>(x) = \J\ogx. □ 


The algebraic part is a consequence of the following result. 

Lemma 14 Let n be squarefree and t be odd. Denote the intersection of Q(Cs) 
and K2nt,2nt by L nt • 

i) If r = 1, 2 { hn, 8|L> and D\8nt, then 


L nt = Q(a/ 2 • sgn(g) ), d(r,tn) = -- c _i(r,tn) = 1 + s S n (^) 


if D = 8(mod 32) and 

L n t = Q(\/2 ■ sgn(-g) ), c x (r,tn) = 1 + s S n (ff) ^ C l ( r tn ) = 1- s S n (s) 


if D = 24(mod 32). 

ii) We have ci(r, tn) ^ C-i(r,tn) if and only if r = 1, 2 { hn, 8|H and D\8nt. 

Proof i). Under the hypothesis of part i) we can apply Lemma 0 to infer that 
L nt , the field intersection in the definition of c±i(l, tn), equals Q(a/ 2 • sgn(g) ) if 
D = 8(mod 32) and Q(y/2 • sgn(— g) ) in the remaining case D = 24(mod 32). 
Writing \/2 = Cs + Cs" 1 one calculates that cr = ti i i(v / 2) = tV%, &±i,i(V—2) = 
i\/ —2. From this observation and the definition of c±\{r,tn) the result follows 
at once. 

ii). ‘<^=’. Follows by part i). ‘=U. The intersection of the fields Q(y+ 2 ) and 
Q(C 2 ^n, y 2rtn ) is abelian and, since 4 { nt, is contained in Q(C 2 H-i, v^2) or 
Q(C 2 h-i, V~2)- As ct i,r (C 2 r + 1 ) = cr -i,r(C2 r + 1 ) for every r > 1, we deduce that 
Qy+z) Fl Q(C 2 r tn, must contain at least one element from {y/—2, V2}. 

Now let us consider how ay ir and a_i ir act on \[2 (= ( 8 + Cs" 1 )- For r > 2 we 
have <7 i i? .(a/±2) = <r_i >r .(A/±2), since then a 1)T .(C 8 ) = cr -i,r(Cs)- Thus we must have 
r = 1. If n is even, then i G Q(y+ 2 ) Fl Q(C 2 r tn, g l ^ rtn ) and since er ±11 (y) = —i, 
we infer that c±\{r, tn) = 0. Thus n is odd. 

From the above discussion it follows that Q(Cs) F Q(C 2 nt,g 1/2nt ) With nt odd 
must contain \[2 or \J—2. By Lemma 0] this leads then to the further restrictions 
(apart from r = 1 and 2 \n)\ 2 \ h, 8\D and D\8nt. □ 


Proof of Theorems |T] and |21 The claims in the cases where the order is even 
or l(mod 2) are already known. The results for the remaining cases are a 
straigthforward consequence of Theorem El and Lemma El If s > 2, then 
r > 2 in the triple sum for the densities. By Lemma El we then have that 
C-i(t, rn) = c\(t,rn) and hence 5 S (1, 2 s ; 1,4) = 5 S (1, 2 S \ 3,4). On noting that 
iV a (l, 2 s ; 1,4)(x) + iV 3 (l, 2 s ; 3,4) = N g (l,2 s ;l,2)(x) the proof of Theorem Q] is 
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then completed. 

Since for ‘most’ t and n we have ci(l,£n) = c_i(l,tn) (Lemma 11411 it is nat¬ 
ural to compute the difference <5^(3,4; 1,4) — S g ( 3,4; 3,4). Since 5 S (3,4;1,4) + 
5 S (3,4;3,4) is easily evaluated, we are then done. We proceed by filling in the 
details. 

Proceeding as in the proof of Theorem El we infer that 


^(3,4; j, 4) 


<5,(1,2; j, 4) 


oo oo 


5 s (l,4;j,4)=X;E 


t= 1 n= 1 

2\t 


h( n ) c y 1 (id( 1 ’ tn) 

[Q(C 8 ,C 2 te ,^/ 2 te ):Q] 


(now only the terms with r = 1 contribute and hence the triple sum is reduced to 
a double sum). From the latter formula and the fact that ci(1, tn) = c_i(l,tn) 
for t is odd and n is even (see Lemma HU), we infer that 


<f g (3,4;l,4)-h 5 (3,4;3,4) = ^^ 

2 ft 2 \n 


hW[cyi(t)(M'n) - c_ Mt) (l,tn)] 
mCs,C 2 tn,g 1/2tn )-.Q] 


On writing nt = v we obtain 


5 fl (3,4;l,4)-5 fl (3,4;3,4) = 5] 

2\v 


V(t>)[ci(l,t<) - C-!(1,d)] 
[Q(C 8 .C2„.9 1/2 *') : Q] 


On invoking Le mm a fTP the latter sum is seen to equal A s /2. We now infer from 
Theorem El that 

N g (3, 4; 1, A)(x) — iV s (3,4; 3,4)(x) = ^Li(x) + O ^log \g!g 2 \ g /2 , 

where the implied constant is absolute. On noting that 

N g { 3,4; l,4)(x) + 7V ff (3,4;3,4)(x) = #{p < x :p = 3(mod 4), (^) = 1}, 

the result easily follows. □ 


5 Explicit evaluation of the densities 


In this section we explicitly evaluate the densities (under GRH), computed in 
Theorem HI and Theorem El 

Using Lemma |U Theorem |l] can be made completely explicit. For reasons of 
space we restrict ourselves to describing the situation for a ‘generic’ g. 


Theorem 7 (GRH). Let s > 2. If h is odd and D contains an odd prime factor, 
then 


<5,(1,2 s ; j, 4) 


2 1-s — | • 4 1-s 

if 3 

4^/6 

if 3 

A 1 - 3 /?, 

if J 


if 3 
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Proof. The conditions on h and D ensure that the degrees : Q] occurring 

in the sums in Theorem Q equal <p(m)n. It then remains to sum some geometric 
series. □ 

Theorem [21 shows that the sets N g (3, 4; j, 4) considered there have a density, 
S g ( 3,4; j, 4), under GRH. The case where j is even is trivial and left to the reader. 


Theorem 8 (GRH). Let g e Q\{ —1, 0,1}, Write g = ±g^, where g 0 > 0 is not 
a power of a rational number. For any prime p define e p by p £p \\h. 

If h is even, then 5 g (3, 4;j, 4) = (1 + sgn(g))/8. 

Next, let h be odd. Then 5 S (3,4;1,4) = <5 g (3,4;3,4) = 1/8, unless D, the dis¬ 
criminant of the quadratic field Q(y/go), is divisible by 8 and has no prime divisor 
congruent to l(mod 4), in which case we have 


1 (-l) 2 ^ 

4( 3 > 4 ;j> 4 ) = q + s g n (^)--o- P 1 P 2 P 3, 


where 


*= n 


p\D 

p=3(mod 4) 


8 


2[p £p - (_!)«*] 2p(—1) 


+ 


pe p -l(^p2 _ 2) p e p(j ) 2 + l)(p — 1) J 


p\D, p\h 
p=3(mod 4) 


2 \p £p - (-l) ep ] 2p(—1) 


P= TT 1 J I 

2 II l 1 T) e P -lfn2_-n 4 


e p +l 


p e p — 1 (j> 2 — l) p £ p(p 2 +1 )(p — i) 
2 p 


, and 


p 3 = n h- 


p\hD 

P =3(mod 4) 


O 2 + i)(p- i)y ’ 


and thus in particular sgn(5 g (3,4; 3,4) — 5 g (3,4; 1,4)) = sgn(g), since the local 
factors of Pi P 2 P$ are all positive. 


Corollary 2 (GRH). For j = 0,1,2,3 we have <5 g (3,4;j,4) = cfij) + C 2 (j)A^ 1 
with Ci(j) > 0 and c 2 (j ) rational numbers. We have c 2 (0) = c 2 (2) = 0 and 
c 2 (l) = —c 2 (3). For a ‘generic’ g all c 2 (j ) will be zero. 

Remark. Using the result going back to Landau that there are 0[xf \Hogx) 

integers n < x having only prime divisors p with p = 3(mod 4) [13, pp. 641-669], 

it is easily inferred that the number of non-generic integers g with \g\ < x is 
0(x/ \/logx). 

Corollary 3 (GRH). Suppose that h — 1. Then S g ( 3,4; 1,4) = <5^(3,4; 3,4) = 
1/8 unless D is divisible by 8 and has no prime divisor congruent to l(mod 4), 
in which case we have 

4(3,4; j, 4) = I + sgn(.g)H Vjl -^-^-JJ —-^ 

8 8 J "B' P ~ P ~ P ~ 1 
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Corollary 4 (GRH). Let j be odd. Then 

^(3,4;j,4) + 5_ ff (3,4;j,4) = i 

An alternative proof of the latter corollary is obtained on combining Propostion 

[T]with Theorem [21 


Proof of Theorem [S] The proof is easily deduced from Theorem [21 

If h is even, then A g = 0 and the result follows on noticing that the density of 
the set of primes p < x with p = 3(mod 4) and (g/p) = 1 equals (1 + sgn(g))/4. 

Next assume that h is odd. Then the set appearing in the formula for 
N g (3, 4; j, 4)(x) given in Theorem [21 has density 1/4. If 8 { D it follows from 
Lemma 01 that the summation conditions are never met and hence A g = 0. So 
we may assume that 8|H. By Lemma[Uthe n’s appearing in the two summations 
are divisible by D/8 and thus if D has a prime divisor p with p = l(mod 4), then 
h^ x (D / 8) = 0 by Lemma [HI and hence h^ 1 (v) = 0 for these v, which shows that 

It remains to deal with the case where 8|D and D contains no prime divisor 
p with p = l(mod 4). Using Lemma El and Lemma [U we infer that 


A 9 = sgn(g)- 


^£±8 
T) 16 


2 tv f h 


/i v , 1 (n)(/r,n) 
v<p(v) 


= sgn (g)- 


^£±2 
T) 16 


C x (h, 2,j). 


On applying Le mm a HIH with r — 2 and s = D/8 (note that since D is a discrim¬ 
inant D/8 must be squarefree), we obtain that 

A i V ^ g +8 P1P2P3 J-|- f ^ f ^ P1P2P3 

A g = sgn(. 9 )(-l) i 6 [[ (- 1 ) = —sgn(g) . 

p\D 

p=3(mod 4) 


An easy analysis shows that the local factors in the products Pi, P 2 and P 3 are 
all non-negative. □ 


6 Modulus 3 

The case d = 3 can be dealt with along the lines of the case d = 4, hence we 
supress most details of the proofs. Our starting point is the following analog of 
Theorem El 

Theorem 9 (GRH). Let £ 0 , £1 be the principal, respectively non-principal charac¬ 
ter modulo 3. For r > 1 let a' lr ,a'_ lr E Gal(Q(^ 3 r+i)/Q) be the automorphisms 
that are uniquely determined by cr l r ((/ 3 r+i) = Q]//)' , respectively cr_ 1 , r (C 3 r +i) = 
C^ :r • For j = —1 and j = 1 let 

d- (r, tn) — l ^ ^ 0 "t,l'l l Q(C3r+l)nQ(C3i'tn;9 1 ^ 3 tn ) ^’ 

1 10 otherwise. 

Let s > 1. For j = 1 and j = 2 we have 

A ff (l,3 s ; j,3)(x) = 6 g (l,3 s ;j,3)-^— + 0 (\og\gig 2 \-—) , 

log a; V log 6/1 xj 
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where 


OO OO OO 




^( n )4i ot)( r ’ tn ) 


r>s £= i n 
“ 3ft 


- [Q^C^n,^*") : 


and the implied constant is absolute. 


It is very easy to see that c\ (r, tn) = c / _ 1 (r, tn) for r > 1 and thus for s > 1 we 
infer that, under GRH, <5^(1, 3 s ; 1,3) = 5 ff (l, 3 s ; 2, 3). Since (unconditionally) 

9 ^ ^ 2 ([AV,3 r : Q] [■^3 r + 1 ,3 r : Q]) 

we then easily deduce, using Lemma [TJ the following result. 


Theorem 10 (GRH). Let e 3 = u 3 (h) and s > 1. If e 3 < s, then 

( 5 g (l, 3 s ; 0, 3) = 3 1 ~ s /2 - 3 2+e3 ~ 2s /8; 

{ 5 g { 1,3 s ; 1,3) = 3 2+e3 - 2 Vl6; 

[ d g (l, 3 s ; 2, 3) = 3 2+e3 ~ 2s /16. 

If e 3 > s, then 

{ 5 S (1, 3 s ; 0,3) = 3 1_e3 /8; 

< <5 fl (l, 3 s ; 1, 3) = 3 l ~ s /4 - 3 1 ” e3 /16; 

{ S g (l, 3 s ; 2, 3) = 3^/4: - 3 1_e3 /16. 

The reason that we cannot take s = 0 in Theorem El is that <j'_ 1 o does not give 
rise to an automorphism of Q(( 3 ). On the other hand a\ 0 does and thus we can 
define c( (0, tn) as in Theorem El 

Let j e {1,2}. Since N g (0,l\j,3)(x) = N g (l,3\j,3)(x) + N g (2, 3; j, 3)(x) + 
0(1), and iV 0 (l, 3; j, 3)(x) is covered in Theorem El it remains to deal with 
N g { 2,3;j,3)(x). Note that 


N g (2,3]j,3){x)= E #{p < x :p = 2(mod 3), r g (p) = t}. 

t=j (mod 3) 

Reasoning as in Theorem El we then hnd that, under GRH, we have 


where 


N g {2,3]j,3)(x) = 5 g (2, 3; j, 3)y—-f O flog \gig 2 \ —) , 

logx V log 6,2 xj 


2,3;j, 3) = E E 


/ n(n)c' 1 (0, tn) 


1 (M* 1 ) +frfo')M t, )) c i(M) 

2 [Q(C3,C»,9 1Al ) :Q] 

i _l ilill sr' Iki ( j? ) 

4 4 ^ [K v v : Q] ’ 

C3&K v ,v 


(13) 


and the implied constant in m is absolute. In the derivation of the second 
equality we used Lemma El and in the derivation of the latter equality we used 
the trivial observation that c\ (0, v) — 1 iff Q 3 ^ K V}V . To sum up we obtained the 
following theorem. 
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Theorem 11 (GRH). The estimate 03 holds with 


S rn o. .• ZoU) , Zl U) 
o 9 (2, 3; j, 3) = —— + —— 


E 

3\v 

C3&K v ,v 


K 0 ) 

[K v>v : Q] 


and an absolute implied constant. 

(Since the condition ^ K v>v implies 3 \ v, the latter condition can be dropped, 
in principle.) Theorem ITT! is the density version of the following lemma. 

Lemma 15 The quantity N g (2, 3; j,3)(x) equals 

<x:p = 2(mod 3)} + ^ hA v )[^K v ,S x ) ~ ^k^,, v {x)]. 

3fu 

Proof. Similar to that of Lemma EH □ 


On noting that 


1 


E 


hA v ) 


2 [ Kv ’ v: 


3fu 


[K v 


[K: 


3v,v 


and invoking Lemma and Lemma 1101 we obtain the following three colloraries 
of Theorem EH 

Corollary 5 (GRH). Put e = 1 if 3 \ D and e = —1 otherwise. If g > 0, then 

n\ 


Wo o. _• o\ Zo(j) , Zi(j) 

s g{ 2,3;j, 3) = —- + y— 

If g < 0 , then S g (2,3;j, 3) equals 


Qi (h, 3,1) + eCfi (h, 3, 


(3, Tli) 


Uj) | Ziti) 




(ni, 3) / 


Corollary 6 (GRH). Recall that e p = v p (h). Define fi(n) = . u pi. n )- Define 


n 


p\D, p >2 
p=2(mod 3) 


2 \p £ p - (-1) 6 ?] 

pe p -l(^p2 _ 1 pe P (^p2 _|_ l)(p _ 1) 


+ 


2 p(-l) e 


and 


n 


p\2D 

p=2(mod 3) 


1 2[p ep - (-l) ep ] 2 p(—l ) ep+1 


p e p~ 1 (p 2 — 1) p e p(p 2 + l)(p — 1) J ' 


Put ei = 0 if D has a prime divisor q that satisfies q = l(mod 3) and e\ = 1 
otherwise. If g > 0, then 

8 g (2,3]j,3) = 3,1) + eh^-{-l) n{ni) 2 e2+2 ~ 2 ^ ni) P[Pf 

If g < 0, then 

„ Uj) , «1 U)c [2« -(-i)«] 

W,3;j,3) = — + — p- 3 . 2 ^-! - 5 -JC £ .(fe,6,1)+ 


Cl 


60') ^ 2e2+2-2u 2 (n 1 ) pi pi 
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Corollary 7 (GRH). Suppose h — 1. We have 


S g (2, 3;j, 3) = + ^l + ei(-l) 0(ni) 2 4 - 2y2(ni) JJ 


2p 


p|D 

p> 3 


p 3 — p 2 — p — IJ 


Remark. Note that 


i) = n (i- x-7r,: + 

p=2(mod 3) 


2[p e * - (-l) ep ] 2p(—1) 


e p +l 


pe p —\ (p 2 _ 1) p e p(p 2 + l)(p — 1) 


A somewhat tedious analysis of Corollary |U] together with Theorem Hill yields the 
following size comparison of 5 S (2, 3;1,3) with <5 S (2,3;2,3) and of 5 ff (l,3) with 
M2,3). 

Proposition 3 (GRH). If g > 0 and h is even, then 8 g (2, 3; 1, 3) < 8 g (2, 3; 2, 3), 
otherwise <5 g (2,3;l,3) > S g ( 2, 3; 2, 3). We have <5 9 (2,3;1,3) = 5 g ( 2, 3; 2, 3) iff 
Q(v / do) = Q(\/3) and, v 2 {h) G {0,2}. The same result holds with S g (2,3;j,3) 
replaced by 8 g (j, 3). 


7 On the generic behaviour of S g (a 1 d ), d — 3,4 

If g is not a square or -1, then an old heuristic model predicts that the number 
of primes p < x such that g is a primitive root mod p should be asymptotically 
equal to ffff, p<x <p{j> — l)/(p — 1), where <p{jp — l)/(p — 1) is the density of primitive 
roots in F*. It is easily proved, see e.g. (20j, that on average ip(p — l)/(p — 1) is 
equal to the Artin constant A, that is 

lim J_y V(P-P =A = Q.37395.... 
ir(x) p — 1 

V 7 p<x 


From the work of Hooley HU it can be deduced that under GRH for a positive 
proportion of all g the above heuristic is false. 

Let 8(p; a, d) = J2 r \ P -i. r=a( mod d) <d( r )/(P “ 1), then S (P', a, d ) is the density of 
elements in the multiplicative group of the finite held F p with order congruent 
to a(mod d). A (naive) heuristic prediction for N g (a,d)(x ) is then provided by 
T < x 8(p; a, d). It can be shown that lim^oo 8{p\ a, d)/7r(x) = 8 (a, d) exists 
fit] . For d = 3,4 some computation j2I| shows that 


{ | if a = 0(mod 3); 

■A + -j 1 - if a = l(mod 3); 
Ye ~ if a = 2(mod 3) 


and 


8(a, 4) 


| if a is even; 
77 if a is odd. 

D 


On comparing this computation with our conditional results for 8 g (a, 3) and 
8 g (a, 4) we obtain the following result. 
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Proposition 4 (GRH). Let d — 3,4 be fixed. There are at most 0(x/\/Togx) 
integers \g\ < x for which 5 g (a, d ) 7 ^ 5 (a, 4) for some integer a. In particular, for 
almost all integers \g\ < x we have S g (a, d) = 5(a, d) for every integer a. 

This proposition shows that for fixed 4 = 3,4 it makes sense to call an integer g 
generic if d g (a, d) = 5(a, d) for every integer a. 

In a similar vein we have: 

Proposition 5 (GRH). Let d = 3,4 be fixed. If \D(g)\ tends to infinity as g 
ranges over a set of rationals g for which h = 1, then 5 g (a, d) tends to S(a, d). 

The latter two results seem to hold for other values of d as well (cf. Table 2 of 
[2TJ ), with Ofx^/\ogx) replaced by o(x). 1 might return to this in a sequel. 
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Vries Institute (Amsterdam) in the NWO Pioneer-project of Prof. E. Opdarn. I 
am grateful to Prof. Opdam for the opportunity to work in his group. 
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(extensive) numerical work. The data given in Tables 1 and 2 were calculated 
using a C ++ program kindly written by Yves G allot. Computations took on 
average (per case) 25 minutes on a 1.50 GHz Pentium IV processor. 
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8 Tables 

We illustrate our results by some examples (assuming GRH). 


Table 1: The case d — 3 


9 

9o 

h 

«l,3)-<5„(2,3) 

numerical 

experimental 

—14 4 

14 

4 

3AJ4, 

+0.13048284... 

+0.13045317 

-196 

14 

2 


+0.17397712... 

+0.17399131 

—3 8 

3 

8 

15^/16 

+0.16310355... 

+0.16310903 

-3 

3 

1 

5A e J2 

+0.43494280... 

+0.43499017 

-2 

2 

1 

37W 8 

+0.06524142... 

+0.06525031 

3 

3 

1 

0 

0 

+0.00001393 

9 

3 

2 

-3AJ2 

-0.43494280... 

-0.43502303 

81 

3 

4 

0 

0 

-0.00001895 

6561 

3 

8 

-bA t Ji 

-0.21747140... 

-0.21748481 

2 

2 

1 

3^,/S 

+0.06524142... 

+0.06515583 

4 

2 

2 

-TA.Ji 

-0.30445996... 

-0.30442279 

5 

5 

1 

67^/94 

+0.12400497... 

+0.12397327 

25 

5 

2 

-15R4 6 /94 

-0.27947388... 

-0.27952119 

49 

7 

2 

-3A t J2 

-0.26096568... 

-0.26097396 

2401 

7 

4 

~A ( J 2 

-0.08698856... 

-0.08697494 


Table 2: The case d = 4 


9 

go 

h 

4,(1,4) -4(3,4) 

numerical 

experimental 

-216 

6 

3 

9A+/28 

+0.20688771... 

+0.20686925 

-9 

3 

2 

0 

0 

+0.00000068 

-81 

3 

4 

0 

0 

-0.00000232 

2 

2 

1 

—Api /4 

-0.16091266... 

-0.16088852 

4 

2 

2 

0 

0 

+ 0.00001122 

8 

2 

3 

—A^, i/28 

-0.02298752... 

-0.02301736 

512 

2 

9 

-3A/J28 

-0.06896257... 

-0.06897632 

216 

6 

3 

— 971^/28 

-0.20688771... 

-0.20687020 

2048 

2 

11 

-489A/J2396 

-0.13136276... 

-0.13134226 

6 9 

6 

9 

/4 

-0.16091266... 

-0.16088478 

6^7 

6 

27 

-23^/84 

-0.17623768... 

-0.17620628 


The number in the column ‘experimental’ arose on taking the density difference 
over the first 10 8 primes, not letting the primes p for which the order of g mod p 
is not defined contribute to either 5 S (1,3) or 5 g ( 2,3) (in Table 1), or 5 S (1,4) and 
<5 S (3,4) in Table 2. Thus, for example, in the column headed ‘experimental’ in 
Table 1 the numbers 

3)(piqs) — N g ( 2, 3)(p 10 s) 

10 8 

are recorded (recall that p 10 s = 2038074743). The last decimals in the columns 
headed ‘numerical’ and ‘experimental’ are not rounded. 
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